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We compute the strong coupling limit of the shear viscosity for the A/" = 4 super- Yang-Mill theory 
with a chemical potential. We use the five-dimensional Reissner-Nordstrom-anti-deSitter black hole, 
so the chemical potential is the one for the _R-charges I7(l)|j. We compute the quasinormal frequen- 
cies of the gravitational and electromagnetic vector perturbations in the background numerically. 
This enables one to explicitly locate the diffusion pole for the shear viscosity. The ratio of the shear 
viscosity r] to the entropy density s is ry/s = l/(47r) within numerical errors, which is the same result 
as the one without chemical potential. 

PACS numbers: ll.25.-w, 11.25.Tq, 11.25.Uv 



I. INTRODUCTION 

Gauge/gravity dualities or AdS/CFT dualities (anti-deSitter/conformal field theory) are particularly interesting at 
finite temperature. One can use the dualities to study the dynamics of gauge theory plasmas or one can use them to 
address the long-standing issues of black holes. 

Finite temperature gauge/gravity dualities were proposed in Refs. 1, 2]. At finite temperature, gauge theories 
are expected to deconfine, so the main evidences of the dualities are rooted in the ones for the deconfinement. It is 
well-known that there are two order parameters for the deconfinement: 

1. Polyakov loops ((P) 7^ in the plasma phase and (P) = in the confining phase) 

2. Large- A^c behaviors of the free energy (P ^ 0{N^) in the plasma phase and F ^ 0(1) in the confining phase) 
As is evident from the nature of these evidences, there are some shortcomings for the finite temperature dualities: 

(i) The dualities are understood only qualitatively. (This is due to the lack of supcrsymmetry.) 

(ii) Thermodynamic properties are compared,^ but it is not evident that the dualities hold even for nonequilibrium 
processes. 

(iii) Gauge theories other than TV = 4 super- Yang-Mills theory (SYM) are less understood. For zero temperature, 
many backgrounds with less supcrsymmetry are known. But even those T = backgrounds are rather involved, 
so their T ^ generalization is a formidable task (See Refs. i2l, iJ, 0, and references therein for some 
exceptions.) 

Remarkably, there is a claim to solve these problems all at once, which is the hydrodynamic description of g auge 
theory plasmas using gauge/gravity dualities Q-liJ- K is originally studied in the context of A/" = 4 SYM 0| and 
later is generalized to many cases. The claim is that gauge theory plasmas which have gravity duals have a universal 
ratio of the shear viscosity rj to the entropy density s at large 't Hooft coupling: 
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This result overcomes the above difBculties in the following sense: (i) The result is a definite quantitative statement; 
(ii) Nonequilibrium processes are considered; (iii) It is expected to be satisfied for a wide range of gauge theories. 
Interestingly, RHIC (Relativistic Heavy Ion Collider) has observed the "elliptic flow" l22| ^ w hich may be a consequence 
of the shear viscosity, and the estimated value is rather close to the duality values ^23, . 

The hydrodynamic description is powerful yet there are many open issues. Of particular interest in view of exper- 
iments is the case with a chemical potential. All discussion so far is limited to zero chemical potential. The general 
theorems of the universality fail in the presence of a chemical potential [^ IllllT^ llSj. However, real experiments are 
done in the finite baryon number density, so one has to consider systems with a chemical potential. 

It is not clear how one can incorporate the baryon number density in the present scheme. However, a charged black 
hole is a simple way to introduce a chemical potential since the electromagnetic potential plays the role of a chemical 
potential for black holes. In this paper, we use the five-dimensional Reissner-Nordstrom-AdS (RN-AdSs) black hole 
to study the issue. 

Our main tool is the numerical computation of the quasinorm al ( QN) frequencies of the black hole. This approach 
has been used to compute the shear viscosity, e.g., in Refs. 0,1^. In hydrodynamics, the shear viscosity accounts 
for the diffusion. In black hole physics, the diffusion is due to the QN modes. 

This is not an easy task however since the perturbation equations are highly-coupled. When the number of degrees 
of freedom are increased, the task is extremely hard. (Thus, many authors often consider the gravitational tensor 
perturbation ^; the tensor perturbation reduces to a minimally-coupled scalar field. Then, one uses the Kubo formula 
to derive the viscosity.) Fortunately, the perturbation equations in question are studied by Kodama and Ishibashi 
|25|: they have decoupled the differential equations, so we use their formalism to obtain the shear viscosity. We 
compute the gravitational vector perturbation and locate the pole explicitly for the shear viscosity. For completeness, 
we also compute the electromagnetic vector perturbation, which is technically very similar to the gravitational vector 
perturbation. 

The plan of the present paper is as follows. First, in the next section, we summarize some background materials; 
We give the relevant physics of the RN-AdSs black hole and discuss the relationship between the QN frequencies 
and hydrodynamics. In Sec. IIIII we present the perturbation equations and a numerical method to obtain the QN 
frequencies. In Sec. II VI our results are discussed. We conclude in Sec. [Vl with discussion and implications. 

After this work is being completed, we learned that Refs. [26L |28| also compute the shear viscosity for charged 
AdS black holes. We briefly discuss their results and compare with our results at the end of this paper. 



II. BACKGROUNDS 

A. Five-dimensional RN-AdS black hole 

The RN-AdSs black hole is a solution of the flve-dimensional Einstein-Maxwell- AdS theory: 
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(2.1) 



(We use the unit SttG = 1 and set the AdS radius I — 1.) The black hole (with noncompact horizon) is given by 

ds^ = -f{r)dt^ + f-Hr)dr^ + r^idxf + dxj + dxj) , (2.2) 
/(r) = r^- 2My + Q^r* = i^' - ^Di^' -J-K^' + ^l) ^ (2.3) 

where Tq = r\ + r^. The black hole with noncompact horizon corresponds to a large black hole limit of the RN- 
AdSs black hole with horizon. 

The surface gravity k and the mass parameter M are written by r_/r^ as 



2M 



(2.4) 
(2.5) 



^ See Sec. Ill Bl for the definition of various perturbation modes. 
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In order to see the gauge theory interpretation of the charge, it is convenient to see the embedding of the black hole 
into ten dimensions. The original finite temperature duality is 

A/" = 4 SYM at finite temperature ^ type IIB string theory in Schwarzschild-AdSs black holes (SAdSs) xS*^ 

One can add an angular momentum along S^, which is known as the "spinning" D3-brane solutions |29l l3Clj| . The 
angular momentum becomes a charge after the reduction. The symmetry of corresponds to an internal symmetry 
of SYM, R-symmetry SO{6). The R-symmetry group 5*0(6) is rank 3, so one can add at most three independent 
charges. The three-charge solution is known as the STU solution [sj. When all charges are equal, the STU solution is 
nothing but the RN-AdSs black hole. In this case, the charge corresponds to the diagonal U(l) of the ?7(l)fj charges. 

Because the charge corresponds to the U{1)r charge, this is by no means realistic. However, the theory has 
interesting features which are common to the real QCD. For instance, the phase diagram 32', 3^] is qualitatively 
similar to the QCD diagram."^ Our system is not a realistic chemical potential, but it may mimic the realistic case 
and one may learn an interesting lesson for the case of chemical potential. 

It is straightforward to compute various physical quantities of the black hole using the Euclidean path integral 
technique j32|]. According to gauge/gravity dualities, they are translated into the corresponding quantities of a gauge 
theory (the energy density e, the pressure p, the temperature Tjj, the entropy density s, the chemical potential ^h, 
and the charge density Q): 

e = 3M, p = M, Th = 7^, s = 27rr^, $h = a/^4' Q = V6Q. (2.6) 

zn V 2 

The speed of sound is then given by Vs = y^dp/de — l/\/3, which results from a scale invariance of the solution (|2.3|l . 
(The symmetry is part of conformal invariance at zero temperature.) The invariance also requires the bulk viscosity 
to vanish C = 0. 



B. QN frequencies and hydrodynamics 

The QN frequencies measure how black hole perturbations decay. (See Ref. [l^ for a comprehensive review.) Since 
the initial value acts as a source, the main interest is the retarded Green's functions with appropriate boundary 
conditions. It turns out that the Green's functions have a number of complex poles in the complex w-plane, which 
means that the perturbations decay exponentially in time. So, the diffusion is mainly governed by the lowest pole. 

To be more specific, let us consider the gravitational perturbations which take the form /ip,y(r) ^-^^t+^kxa _ -y^g look 
for the lowest pole which survives in the low-frequency, long- wavelength limit u^k) ^ 0, fc — > (hydrodynamic limit). 
The diffusion is due to the black hole absorption. One interprets the diffusion as a consequence of viscosities of the 
dual gauge theory plasmas. 

The perturbations can be decomposed by the world- volume symmetry 0(3) of xi , • • • , ^3. (We follow the conventions 
of Ref. H^.'*) Then, the gravitational perturbations are decomposed as the tensor mode, the vector mode ("shear 
mode"), and the scalar mode ("sound mode").^ 

Such a decomposition is essentially the same as hydrodynamics. In hydrodynamics, the vector mode measures the 
shear viscosity, and the scalar mode measures the bulk viscosity and the speed of sound. For example, as we will see 
below, the vector mode satisfies the following dispersion law (in the hydrodynamic limit): 

^^-.dJ'-)\ A,:^^. (2.7) 

Then, one can regard 77 as the shear viscosity of the dual gauge theory plasma. 

Similarly, the electromagnetic perturbations are decomposed as the vector mode and the scalar mode. In this paper, 
we compute the QN frequencies for the vector modes since we are mainly interested in the shear viscosity. 

As is evident from the above discussion, this is a Lorentzian computation. On the other hand, gauge/gravity 
dualities are originally defined in the Euclidean signature (such as the GKP-Witten relation P,0]). Many authors have 



^ This is the case of the RN-AdS black hole with compact horizon or compact SYM. This paper considers the RN-AdS black hole with 
noncompact horizon, which is always in plasma phase. 

This convention is slightly different from the one of Ref. |21| . In Ref. |21| . the perturbations /i^i,(r)e~'"'+*''^3 are considered and they 
are decomposed by the remaining world- volume symmetry 0(2) of xi and X2. 
^ In Ref. I21I . these modes are called the spin 2, spin 1, spin modes, respectively. 
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studied the Lorentzian prescription llSllSSllSft 37] and we essentially use the one given in Ref. fl5|. In principle, one 
can use only the Euclidean prescription and analytically continue the results to the Lorentzian signature. However, 
the differential equations in gravity backgrounds often have more than three regular singularities, so the analytic 
computations are not possible; one can locate only the pole locations numerically. Without knowing the full analytic 
properties, the analytic continuation is not possible. So, one has to work with the Lorentzian signature. In simple 
cases where the analytic solutions are available (such as the three-dimensional black hole) , the Euclidean prescription 
and the Lorentzian prescription by Ref. 113] give the same results. 



III. NUMERICAL METHOD 

The perturbative equations for the Einstein-Maxwell- AdS system H2.1|l have been studied in Ref. (See Sec. 4.3.1 
of the paper for the vector perturbations, which are our primary interests.) Let us outline their procedure briefly. 
One first has to find the gauge-invariant variables from the gravitational and electromagnetic perturbations. Denote 
the gauge-invariant gravitational and electromagnetic vector perturbations expanded in vector harmonics by and 
A., respectively. The equations of O and A themselves are not very useful however since they are highly coupled. But 
take the following linear combinations of the gauge-invariant variables: ^ 

^±^a±r'^^^n + b±r^/^A. (3.1) 

Then, the differential equations decouple from each other. In the limit of zero charge, the fields $- and $+ reduce 
to the pure gravitational mode and pure electromagnetic mode, respectively. 

Let us use the dimensionless coordinates f = r/r+ and f_ = r_/r+. The decoupled equations are then given by 

/ W 4 4 + - V{r)] $± (f ) , (3.2) 

(3.3) 



dr df 

3 . 39 Q2 

■ r' ' 
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fi±:=-—M±A, A := V64M2+T2fc2Q2 , (3.4) 

where 2M = 1 + ft + rt and = + 

The QN modes are the solutions of the perturbative equations with appropriate boundary conditions. The boundary 
condition at the outer-horizon is the purely ingoing mode, which is physically natural for the retarded Green's 
functions. The boundary condition at infinity is chosen to be a decaying mode. The asymptotic solutions behave as 
$ ^ f~3/2 Qj. (j) ^ and the former one is chosen by regularity. Incorporating the boundary conditions both for 
the horizon and infinity, we obtain the form 

<f±(r) = (3.5) 



The equation ()3.2|l is transformed into the form 



= ( s{z) t{z) Z—+ u{z) ] iP±{z) , (3.6) 



^ The system studied by Kodama and Ishibashi is the general Einstein-Maxwell system with a cosmological constant in d = n + 2 
dimensions. Set n = 3, K = (the spatial curvature of the horizon), and A = — 1 (the sign of the cosmological constant) in their 
equations. Also, see the paper for the definition of the numerical coefficients a± and b±. 
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where z = (r^ — — r^) and 

s{z) = 4(-l + fl){2 + + z){-l + z){ - 2 + z + f2 (-1 + 2z) Y 

=4(2 + f2)| _2 + z + f2 (-l + 2z) } 

(-1 + + f2 ) ( - 2 - 4(-2 + z)z + 2flz(-l + z + z^) 



l + z{ 3 + (-7 + z)z } ] ^ 



+ ic^(-l + z)(-l + f^z){ -2 + z + f^(_i + 2z) } 



(3.7) 



(3.8) 



u[z) — z 



- (-1 + f2 )(2 + f2 )2| _ 2 + z + f2 (-1 + 2z) } 



X fc2(-l + f'iz) + (-1 + z){ 4(1 + f2 + f4 _^ ^ I ^ 
+ 2iw(-2 + f2 + f4)| _2 + z + f^(-l + 2z) }{ -5 + 3z + f2 (_3 + 52) } 
+ (-l + f^z)(^ 8+(-5 + z)z + f^z2 + f'^z(-3 + 5z) + 3fi{ 1 + z(-5 + 3z) } 



f2 



2 { 10 + z(-19 + 5z) } 
The solutions are expanded as a series to solve the differential equation: 



OO 

n=0 



Then, the expansion coefficients a„ satisfy the following recurrence relations: 

aofli + /3oao = 0, 

aia2 + /3iai + 7109 = 0, 

Oi2az + /32a2 + 7201 + 62aQ = 0, 

anO-n+l + 0nan + lna,n~l + <5na„-2 + PnOn-S =0, n = 3, 4, 



(3.9) 



(3.10) 



(3.11) 



where the recurrence coefficients, which depend on w, are given in appendix. The purely ingoing wave boundary 
condition at the horizon is reflected by the first equation in H3.11|l . and the fall-off condition <I> ~ f-3/2 infinity 
is satisfied for the series il]±{z) which converges uniformly as z — > 1 (namely, for the series where X) '^n exists and is 
finite.) Both conditions are satisfied by special complex values of lo which are the QN frequencies. 

The convergence problem of the series translates into the convergence of the continued fraction (in a compact 
notation of continued fractions) 



ai 



-ii ai72 "273 



Here, {a^}, {/3^}, and {7^} (n = 0, 1, 2, • • • ) are the coefficients of the three term recurrence relation 

a'^ai + /3oao = 0, 

a'^an+i + /3^a„ + 7^a„_i = 0, n = 1, 2, • • • , 



(3.12) 



(3.13) 



which is obtained from the original recurrence relation (|3.11|) by the Gaussian elimination. Thus, the QN frequencies 
are obtained by solving the equation 38] 



^ ^ -il «l72 «273 
Oc'o P'l- P'2- 03- ■■ 



(3.14) 
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IV. RESULTS 

We compute the QN frequencies for <i>±, which are hnear combinations of the gravitational and electromagnetic 
vector perturbations. We first discuss general structure of these poles and thier gauge theory interpretations. Then, 
we closely study the pole which survives in the hydrodynamic limit. 

A. General QN spectrum 




-6-4 -2 2 4 6 -6 -4 -2 2 4 6 

Re((B/K) Re(K)/K) 



FIG. 1: (color online). QN spectrum of the perturbations $_ (left panel) and "1>-|- (right). Dependence of uj/2-KTYi on 
r-/r+ = 0(x), 0.5(+), and 0.95 (*) is shown for spatial momentum k/2nTH = 1.0. On the left panel, the poles which are 
enclosed by the elongated circle are the hydrodynamic pole discussed in Sec. liV 51 In contrast, there is no such poles for $+. 
All poles stay at a finite distance from the real axis. 

The poles for the perturbations $± are shown in Fig. (<i>_ on the left panel and on the right panel). On the 
left panel, the poles which are enclosed by the elongated circle are the hydrodynamic pole; they approach the real 
axis in the limit uj{k) ^ 0, k 0. We closely study the pole in Sec. lIVBl In contrast, ^+ does not have a pole which 
survives in the hydrodynamic limit. All poles stay at a finite distance from the real axis even in the limit. 

Let us look at gauge theory interpretation briefly. According to the dictionary of gauge/gravity dualities, the 
gravitational field couples to the energy-momentum tensor of the gauge theory, so the gravitational perturbations 
correspond to the energy-momentum tensor correlators. The electromagnetic field couples to the i?-current, so the 
electromagnetic perturbations correspond to i?-current correlators in the dual gauge theory side. 

In our choice of field variables the gravitational and electromagnetic perturbations are entangled with each other 
[See Eq. H3.1|l ]. So, our computation mixes the energy-momentum tensor correlator with the i?-current correlator. In 
general, one has to go back to the original variables to extract information of gauge theory correlators. However, in 
this case, such a disentanglement is not really necessary in order to know the plasma viscosity. This is because has 
no hydrodynamic pole. The dispersion relation of the shear mode is simply determined by the hydrodynamic pole 
of This also means that the vector mode of the i?-current correlator has a hydrodynamic pole with the same 
dispersion relation as the energy-momentum tensor correlator. [Namely, there is a pole which satisfies the dispersion 
relation H2.7|l with the same diffusion constant as D^.] 

The poles beyond the hydrodynamic pole are not relevant to hydrodynamics but are interesting from the point 
of view of gauge/gravity dualities. One feature of these poles is that they lie approximately linearly. The slope is 
getting sharper as r_ is increased. This property that poles lie linearly has been observed for the SAdS black holes 
(see, e.g., Refs. [HlHlal^) and for the RN-AdS4 black hole This is also observed for the near-horizon limit 
of nonextreme Dp-branes |42(|. However, its origin and significance are far from obvious. 
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r_/r+ 


a 


b 


c X 10" 




0.0 


-0.500051 


2.00002 


-3.08339 


1.000102 


0.05 


-0.498175 


2.00002 


-2.7581 


1.00010035 


0.1 


-0.49255 


2.00002 


-2.94714 


1.00010002 


0.15 


-0.483179 


2.00002 


-2.51994 


1.00009449 


0.2 


-0.470089 


2.00002 


-2.54955 


1.0000913 


0.25 


-0.453336 


2.00002 


-2.33166 


1.00008669 


0.3 


-0.433032 


2.00001 


-1.86945 


1.00007874 


0.35 


-0.409367 


2.00001 


-1.62585 


1.0000732 


0.4 


-0.382616 


2.00001 


-1.39078 


1.0000651 


0.45 


-0.353155 


2.00001 


-1.17109 


1.00006151 


0.5 


-0.321446 


2.00001 


-0.971175 


1.00005422 


0.55 


-0.288032 


2.00001 


-0.791317 


1.00004974 


0.6 


-0.253502 


2.00001 


-0.628083 


1.00004391 


0.65 


-0.218464 


2.00001 


-0.478805 


1.00003943 


0.7 


-0.183507 


2.00001 


-0.343142 


1.00003295 


0.75 


—0 149179 




—0 293444 


1 00009414 


0.8 


-0.115927 


2.0 


-0.126007 


1.00001675 


0.85 


-0.0841502 


2.0 


-0.0553647 


1.00001171 


0.9 


-0.0541242 


2.0 


-0.0130611 


1.00000517 


0.95 


-0.0260391 


2.0 


0.0549852 


0.999999434 



TABLE I: Fitting parameters and 7num. for various charge-to-mass ratios. With a good accuracy, b = 2, c = 0, and 7n 



B. Hydrodynamic pole 



Let us look at the hydrodynamic mode of in more detaiL In the hydrodynamic Hmit uj(k) ^ 0, A; ^ 0, the 
gravitational vector mode (shear mode) is expected to satisfy the following dispersion law: 



- = -iD„ - 



A, := 



K ' \K J ' e+ p 

We numerically estimate the diffusion constant Z),, through the QN frequencies. The ratio rj/s can be written as 

7] €+p 



(4.1) 



K S 

2(e + p) 



1 



47r 
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„ _ 4[l + (r_/r+)^ + (r_/r+)4] 

Ths " 2-(r_/r+)2-(r_/r+)4 



(4.2) 
(4.3) 



where Eq. 1)2. 6|l is used. The quantity j is rj/s in imits of 1/ (47r). We estimate 7 below; 7=1 means that rj/s is the 
same as the one without the chemical potential. 

The diffusion constant D,, is estimated by choosing the following fitting function: 



— = a X - + c , 



(4.4) 



where a, h and c are the fitting parameters. Then, the parameter a is the numerically estimated value of Z?^ and the 
numerically calculated 7 is given by 



In 



4 [l + (r_/r+)2 + (r_/r+)^ 



(4.5) 



2-(r_/r+)2-(r_/r+)4 

The result is summarized in Table 13 and Figs. El El As one can see from the last column of the table, 7num. = 1 
with a good accuracy, which means that 

1] 1 
s 47r 



(4.6) 
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FIG. 2: (color online). Dependence of cj//27rrH on {k/2'KTYif for r_/r+ = 0(x), 0.25(+), 0.5(*), 0.75(Q), and 0.95(B). The 
slope determines the diffusion constant D,, . 
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FIG. 3: (color online). Dependence of four low-lying QNM's on k/2nTH for r-/r+ — 0(x), 0.5(+). The left panel shows the 
real part of uj/2nTu and the right shows the imaginary part. The hydrodynamic pole is the one which approaches — > as 
fc — > (the bottom one on the left panel and the top one on the right panel). As one can see, the other poles stay at a finite 
distance from the origin. 



even in the presence of a chemical potential. This is our main result. 



V. DISCUSSION 



We found that rj/s is again l/(47r) in the presence of a chemical potential. A nonzero value of viscosity is probably 
not very surprising for fluids at finite temperature. There is an intuitive argument why the viscosity bound exists 
|12|. Generically, the viscosity is given by 77 ^ eTmft, where e is the energy density and r^ft is the mean free time. The 
entropy density is roughly the order of the number density n. So, 77/s is the order of average energy of the particle 
times the mean free time. From the uncertainty relation, this product must be larger than h. 

This argument is an intuitive argument, but this suggests that nonzero viscosity is also expected in the presence 
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of a chemical potential (although one cannot learn the precise value from the argument). Our computation confirms 
the naive expectation and 77/s is the same as the one without a chemical potential. 

Such a dimensional argument should also hold for the other transport coefficients. So, one would expect some 
relations among various coefficients. (In some cases, one would expect a bound similar to the viscosity bound.) In 
fact, the existence of the relations among various coefficients is well-known in the context of the linear response theory. 

In thermodynamic description, the viscosity bound also constrains the behavior of the relaxation time Trciax- From 
V Cmft and e ^ nT, one obtains rj/s Tr^ff So, the bound implies that Tmft is independent of the coupling in 
the strong coupling limit and depends only on the temperature. Since Tioiax > ''mft, the saturation of Tmft mean that 
Troiax cannot be small indefinitely either. 

In fact, Ref. points out that the relaxation time is saturated at strong coupling and is the order of thermal 
wavelength. The relaxation time becomes shorter as the coupling is increased, but cannot be smaller than the time 
scale determined by the thermal wavelength and causality. In gravity description, this is a consequence of a scale 
invariance (which is part of conformal invariance at zero temperature) for the large black hole limit of SAdS black 
holes. The noncompact {p+ l)-dimensional SYM, which is described by the near- horizon limit of Dp-branes in gravity 
side, are no longer conformal, but this is also true for these systems [42 |. 

The hydrodynamic description using gauge/gravity duals has many open issues. For example, there are several 
independent methods to calculate the viscosity such as the Kubo formula arguments, "membrane paradigm" -like 
approach, and the direct computation of the QN frequencies. Currently, these three gave the same results (for the 
cases where the comparison is possible), which justify these approaches, but the relationship is not clear. Basically, 
these methods can be classified into two groups: 

1. Phenomenologically, the viscosity is given by the one-point correlator of the energy- momentum tensor. 

2. Microscopically, the viscosity is given by the two-point correlators of the energy-momentum tensor, namely by 
the Kubo formula. 

These two methods are of course not independent. Rather, they are related by a linear response theory. Thus, one 
may be able to establish the connection among these methods if one really establishes the linear response theory for 
AdS black holes. This in turn may give the nonequilibrium verification of finite-temperature gauge/gravity dualities. 

After this work is being completed, we learned that three other groups also compute the shear viscosity for charged 
AdS black holes 26, 27, 28]. For the convenience of the readers, we summarize their results and the differences 
from our paper. All four groups have obtained rj/s — l/(An) even with chemical potentials, but the systems and the 
techniques are slightly different from each other: 

• Ref. 26] studied the five-dimensional STU black hole. The gravitational tensor perturbation is computed 
analytically and the shear viscosity is determined by the Kubo formula. 

• Ref. |27] studied the five-dimensional STU black hole. The gravitational tensor perturbation is computed 
analytically and the shear viscosity is determined by the Kubo formula. For the STU solution with a single 
charge, the lowest QN mode of the gravitational vector perturbation is also computed analytically. 

• Ref. ^2§\ studied the RN-AdS4 black hole. The lowest QN mode of the gravitational vector perturbation is 
computed numerically (for a small chemical potential). The shear viscosity is expanded as a power series in the 
chemical potential, and no correction was found for the first few series. 

• We studied the RN-AdSs black hole. The QN modes of the gravitational and electromagnetic vector perturba- 
tions are computed numerically. 
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APPENDIX: RECURRENCE COEFFICIENTS 



The recurrence coefficients in Eq. (|3.11|) are given by 
a„ = -4(1 + n){2 + r'if{K + nK -ioj) , 



13,, ^ (2 + rl) 



K(2 + f2 ) 

12 + /fc2 + An{ 5 + 4n + (-2 + 7n)f^ + (_3 + „)f4 } ± A - 4k 
2iuj{ 2n{2 + rl){4: + 7P'i+ri) + nib + Srl) } - lu'^ {4 + 3rl) 



(A.l) 



(A.2) 



7„ = -K (2 + f2 ) 1^ P(l + 4f2 + fi) + 4 3 + n + + { 16 + n{-29 + 22n) )fi 
+ (-1 + n)(-37 + 22n)fi + (-1 + ?i)(-16 + hn)rt 
+ 3(1 + f^) ± A - 12(1 + f-_)K ) 



2i LO 



K (11 + 26fl + IIP ) + 2(-l + n)(l + fl)(2 + fl)(5 + 17fl + 5P ) 



+ cj^ (5 + 27fi + 25P + 6P ) , 

(5„ = K (2 + + 2fi) [ 4 + (96 + fc^)f^ ± A - 4k 

+ 4{ (-1 + n)n + n(-26 + 7n)f^ + (-2 + n)(-13 + 4n)fi } 
-2iw(l + 2f2)| 2(-2 + n)(2 + f2)(i + 7f2 _|_4^4-j_|_(-3_j_5^2-j^ | 

- (1 + 4f^ + 2fi)(l + 6f?. + 2ft) , 



4k (-3 + n)2(2 + f2 )(i + )2 

+ 4iw(-3 + 7l)(2 + f?.)(l + 2f^)2 +^2 (^-^ _^ _^ g^4 _^ _^ ~8 ^ 



(A.3) 



(A.4) 



(A.5) 
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